Astronomy & Astrophysics manuscript no. mancini 


February 2, 2008 


(DOI: will be inserted by hand later) 





LMC Self-lensing from a new perspective 

L. Mancini^'^'"^, S. Calchi Novati\ Ph. Jetzer\ G. Scarpetta^'^'^ 



o 
o 

(N 



00 

(N 
> 

in 
o 

o 



Oh 
I 

o 



X 



^ Institut fiir Theoretische Physik der Universitat Zurich, CH-8057 Zurich, Switzerland 

^ Dipartimento di Fisica "E.R. Caianiello", Universita di Salerno, 1-84081 Baronissi (SA), Italy 

International Institute for Advanced Scientific Studies, Vietri sul Mare (SA), Italy 
* Istituto Nazionale di Fisica Nucleare, sez. Napoli, Italy 

Received / Accepted 

Abstract. We present a new analysis on the issue of the location of the observed microlensing events in direction 
of the Large Magellanic Cloud (LMC). This is carried out starting from a recently drawn coherent picture of the 
geometrical structure and dynamics of the LMC disk and by considering different configurations for the LMC bar. 
In this framework it clearly emerges that the spatial distribution of the events observed so far shows a near-far 
asymmetry. This turns out to be compatible with the optical depth calculated for the LMC halo objects. In 
this perspective, our main conclusion, supported by a statistical analysis on the outcome of an evaluation of the 
microlensing rate, is that self leasing can not account for all the observed events. Finally we propose a general 
inequality to calculate quickly an upper limit to the optical depth along a line of view through the LMC center. 
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1. Introduction 

The microlens ing; surveys towards the Large M agellanic 

Cloud (LMC) l)Alcock et ani2000aHLasserre et al 200Q) 

have demonstrated the existence of compact objects that 
act as gravitational lenses somewhere between us and the 
LMC. In some cases the distance and the mass of the 
lenses have been determined, thanks to the proper motion 
of the lens observed by the Hubble Space Telesco pe (HST) 
l)Alcock et al. ll2001atlGould ..2004tlDrake et al. ..2004.') or 
to the additional informatio n carried by binary systems 
llAlcock et al. Il2000h[ EooTbh . However, these are special 
cases since for most events only the duration and the posi- 
tion on the sky plane have been measured. These informa- 
tion are not enough to establish definitively if the detected 
events are really caused by white dwarfs or MACHOs in 
the halo of the Milky Way (MW), or are due to stars or 
MACHOs within the LMC itself. 

The survey of the MACHO team indicates a most 
probable Galactic halo fraction of 20%, with limits of 5% 
to 50% at the 95% confidence level, assuming that all the 
events are due to halo lenses. The preferred value for the 
lens mass is ~ 0.4 Mq. This is consistent with the EROS 
survey results, that are given however as an upper limit 
for the Galactic halo fraction. 

An interesting alternative is that of "self lensing" , 
where both source and len s belong to t he lu minous par t 
of the LMC as suggested bv lSahu 1 111991 and lWu I lll994l) . 
However, the initial estimates of the optical depth and mi- 



crolensing rat e were lower than the measured one llGould I 
Il995; Alcock et al. Ill997l EoOOab The self-lensing expla- 
nation has been further analyzed going beyond the hy- 
pothesis of a "simple" geometry for the LMC with disk 
and bar coplanar, so that their relative distance would 
enhance t he optical depth and, t herefore, the rate. In the 
model of IZhao fc Evans ] (pOOff) the disk and bar stars 
are on two distinct planes with different inclinations, so 
that stars on the front plane could lens those in the plane 
« 1 kpc behind. Besides the morphology, another aspect 
considered has been the dynamics of the l umino us compo- 
nents within the LMC. lAubourg et al. I l|l999j) . by using 
a model which takes into account the correlation between 
the mass of the lenses in the LMC and their velocity dis- 
persion, have been able to reproduce a self-lensing optical 
depth, event rate and event duration distribution com- 
patible with the observed ones. Yet, objec tions to this 
mode l have been raised by d ifferent authors ijGvuk et al. I 
I2OOOI: lAlves fc Nelson I l2000l) . especially with respect to 
the adopted distribution and velocity dispersion of the 
lensing stars, which seem to be inconsistent with the ob- 
servations. 

The analysis of Jetzer et al., (2002, hereafter Paper 
I) has shown that probably the observed events are dis- 
tributed among different components (disk, spheroid and 
galactic halo, the LMC halo and self-lensing). This means 
that the lenses do not belong all to the same population 
and their astrophysical features can differ deeply from one 
another. 
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In this paper we address once more the question of 
the presence of a self-lensing component within the LMC 
itself. To this end a correct knowledge of the structure 
and dynamics of the luminous components (disk and bar) 
of the LMC is essential. Here we take advantage of some 
recent studies of the LMC disk (see Sect. 12 .ill , while we al- 
low for different geometries for the still poorly known bar 
component, to calculate the main microlensing quantities. 
Moreove r, with respect to Pap er I, based on the moment 
method I de Ruiula et al. Il99ll) . we perform instead a sta- 
tistical analysis starting from the differential rate of the 
microlensing events. 

The paper is organized as follows: in Sects. |2 and |3| we 
discuss the LMC morphology and present the models we 
use to describe the spatial density of the MW halo and of 
the various components of the LMC. Sect. ^ is devoted to 
the calculation of the microlensing quantities, the optical 
depth and the microlensing rate, as well as to a statistical 
analysis of the self-lensing events. In Sect. El we discuss 
the spatial asymmetry with respect to the line of nodes of 
the observed microlensing events. An improved inequality 
for the optical depth for self lensing by a stellar disk is 
derived in Sect.El We conclude in Sect. [7| with a summary 
of our results. 

2. Observational data 

2.1. LMC Disk morphology 

Recently, in a series of thr ee interesting papers 

llvan der Marel &: Cioni I l200lt Ivan der Marell l200ll 
Ivan der Marel et al. 1 120021) . a new coherent picture of 
the geometrical structure and dynamics of the LMC 
disk has been given. By usin g the data frorn two 
near-infr a.red surveys ^ENIS ( Cioni et al. 1 120001) and 
2MASS l)Cutri et al. I mOO.) stellar catalogs, the LMC 
geometry has been constrained by carefully studying 
the star count map and the characteristic apparent 
magnitudes of Asymptotic Giant Branch (AGB) and 
Red Giant Branches (RGB) stars. The results have been 
further improved by a detailed and sophisticated study 
of the LMC line-of-sight velocity field, coupled with a 
multi-parameter fit on the available velocities of 1041 
carbon stars. A basic assumption is that the carbon star 
population is representative for the bulk of the LMC 
stars. 

The first important conclusion is that the intrinsic 
shape of the LMC disk is not circular, as assumed before, 
but elliptical, with an intrinsic ellipticity e" = 0.312 ± 
0.007. The inclination angle of the LMC disk plane is 
i = 34.7° ± 6° and the hne-of-nodes position angle is 
e = 129.9°±6°. This value is quite different from the LMC 
disk major axis position angle, PAmaj = 189.3° ± 1.4°, 
corresponding to a position angle $ = 202.7° ±1.9° when 
measured in the equatorial plane of the LMC disk, start- 
ing from the axis pointing towards the North. The radial 

^ Near-infrared surveys provide a very accurate view, thanks 
to their insensitivity to the dust absorption. 



number density profile along the major axis follows, to 
lowest order, an exponential profile with an intrinsic scale 
length equal to 1.54 kpc. 

A second important conclusion is that the center of 
mass (CM) of the carbon stars is consistent with the cen- 
ter of the bar and with the center of the outer isophotes of 
the LMC. As a consequence, the idea of using the distribu- 
tion of neutral gas as good tracer for the disk stars, that 
leads to an incorrect LMC model, must be abandoned. 
The obtained values of the right ascension and declination 
of the CM, given in J2000, are acu = 5^ 27.6™ ± 3.9™, 
ScM = -69.87° ± 0.41°. The weighted mean of the ro- 
tation velocity in the range 4-8.9 kpc, where the rotation 
curve is approximately flat, is y = 49. 8± 15.9 km/s, about 
40% lower than the previously inferred and accepted value. 
Taking into account the asymmetric drift effect, the circu- 
lar velocity has been corrected and estimated to be equal 
to V^irc — 64.8 ± 15.9 kms'^. The line-of-sight velocity 
dispersion has an average value a = 20.2 ± 0.5 km s"'^, 
and is little dependent on the radius. The rate of inclina- 
tion change is = — 0.37±0.22 masyr"^, a value similar 
to that determined from N-body simulations by Weinberg 
(2000), which predicts the LMC disk precession and nuta- 
tion to be due to tidal torques generated by our Galaxy. 

A third important conclusion is that the LMC disk 
has a considerable vertical thick ness, in agreemen t also 
with the numerical simulations of I Weinberg I l|200n() . The 
thickening of the LMC disk is due to the gravitational 
interaction with the MW. The ratio = 2.9 ± 0.9 is even 

cr 

lower than the corresponding value for the MW thick disk 
(■^ ~ 3.9). The LMC disk is also flared. The best fit of 
the observed velocity dispersion profile with isothermal 
disk models, whose vertical density profi le is proportional 
to sech^(^), confirms the result found bv lAlves fc Nelsonl 
(j200Q} that the scale height must increase with radius. The 
vertical thicken i ng is also in agreement with the results of 
lOlsen fc Sal^(|2002l) . who argued that the LMC contains 
features that extend up to 2.5 kpc out of the plane. 

Let us note that rece ntly some of these con clusions 
have been challenged by iNikolaev et al. I l)2004l) . whose 
analysis is based on a combination of the results of the 
MACHO collaboration on the LMC Cepheides with the 
2MASS All-Sky Release Catalog. 

2.2. LMC Dark Halo 

Another important conclusion of the van der Marel et al. 
analysis, is connected with the LMC dark halo. Comparing 
the LMC total mass enclosed in a sphere of 8.9 kpc, dy- 
namically inferred to be (8.7 ± 4.3) x 10^ Mq, with re- 



2.7 X lO^Mrr^, and to th e 
" Kim et a,1. 1(1993) . 



spect to the visible mass M^i 

mass of neutral gas « 0.5 x IO^Mq ^ ^, 

Ivan der Marel et al. 1 1)2002|) have concluded that the LMC 
must be surrounded by a dark halo with a radius equal 
to the LMC tidal radius rt = 15.0 z t 4.5 kpc, a valu e 
higher than that, 11 kpc, estimated bv lWeinberg I l|200nl) . 
Moreover, at a distance of 5 kpc from the center of the 
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LMC, the ratio between the tidal force and the LMC self 
gravitational attraction is reduced by 20%. Thus, one ex- 
pects that tidal effects influence deeply the structure of 
the LMC halo, causing probably an elongation of it. 

3. Models 

We adopt the same coordinat e systems and notations as 
in Ivan der Marel et al. I l)2002f) . The position of any point 
in space is uniquely determined by the two angles, right 
ascension a and declination 5, and by its distance D. We 
introduce also a cartesian coordinate system (x, y, z) that 
has its origin O in the center of the LMC, whose coordi- 
nates are (ao, (^Oj -Dq), with the x-axis antiparallel to the 
right ascension axis, the y-axis parallel to the declination 
axis, and the z-axis pointing towards the observer. To 
describe the internal kinematics of the galaxy it is oppor- 
tune to introduce a second cartesian coordinate system 
(x', 2/', z') that is obtained from the system (x, y, z) by a 
counterclockwise rotation around the z-axis by an angle 
equal to the position angle of the line of nodes — 
followed by a rotation around the new a;'-axis by an angle 
equal to the inclination angle of the disk plane i. With this 
definition the (x', y') plane coincides with the equatorial 
plane of the LMC galaxy, and the x'-axis is along the line 
of nodes. 

3.1. LMC Disk 

Following the new results on the LMC morphology dis- 
cussed in Sect. 121 we consider an elliptical isothermal flared 
disk tipped by an angle i = 34.7°±6.2° with respect to the 
sky plane, with the closest part in the north-east side. The 
center of the disk coincides with the center of the bar lo- 
cated at J2000 (a, S) = (5'*27.6" ± 3.9"^, -69.87° ± 0.41°) 
and its distance from us is Dq = 50.1 ± 2.5 kpc . Here 
we take a bar m ass Mbar — l/5Mdisk l|Sahu I Il994l: 
iGvuk et al. Il200(1) . with the visible mass in disk and bar 
Mhar + Mrii.k = (2.7±0.6)x lO^M© l|van der Marel et al. I 

The vertical distribution of stars in an isothermal disk 
is described by the sech^ function, therefore the spatial 
density of the disk is modelled by: 



Pd 



sech 



c 



Cd(i?) 



(1) 



where q = 0.688 is the ellipticity factor, i?d — 1.54 kpc is 
the scale length of the exponential disk, R is the radial 
distance from the center on the disk plane. Cd{R) is the 
flaring scale height, which rises from 0.27 kpc to 1.5 kpc at 
a dist ance of 5.5 kpc from the center ( van der Marel et al. I 
I2OO2I) . and is given by 



Cd{R) = 0.27+ 1.40 tanh 



Af = 0.2765 is a normalization factor that takes into ac- 
count the flaring scale height. 



The cartesian coordinates (^, 77, C, = z') are obtained 
from the system (a;', y' , z') by rotation around the z'-axis 
by an angle equal to (<& — — .|), where $ is the position 
angle of the LMC disk major axis. In this way the £, — rj 
plane coincides with the equatorial disk plane and the r] 
axis is directed along the major axis of the elliptic disk. 

The velocity dispersion is a crucial parameter for the 
estimate of the microlensing rate. The kinematic studies 
of the LMC disk have shown that measurements of the 
velocity dispersion along the line of sight vary between 
roughly 6 and 30 km/s, according to the age of the tracers 
and show little var iation as a function of the radius (see 
iGviik et al. I l|2000tl where a comprehensive table is given). 
In particular, younger populations have a smaller velocity 
dispersion than older ones, a s in the MW. In the picture of 
Ivan der Marel et al. I l(2002l ) within a distance of about 3 
kpc from the center of the LMC, the line of sight velocity 
dispersion (evaluated for carbon stars) can be considered 
as constant with crd,ios — 20.2 ±0.5 km/s. This represents 
our choice for the line of sight velocity dispersion of the 
disk stars. 



3.2. LMC Bar 

An optical bar of roughly 3° x 1° angular size stands out 
in the central region of the LMC. Unlike the disk, the 
geometry of the bar region, in particular its vertical struc- 
ture, is not well constrained. Understanding how the bar 
is displaced with respect to the disk is one of the last 
information needed to complete the puzzle of the LMC 
morphology. 

In Paper I we have described the bar by a gaussian 
density profile following iGviik et al. I ll200nl) . In this paper 
we choose, instead, a bar sp atial density that take s into 
account its boxy shape, as in lZhao fc Evans 1 1|2000|) : 



2Mb 



Ph 



' Rl Sb 



(2) 



where Sb — 1.2 kpc is the scale length of the bar axis, 
i?b — 0.44 kpc is the scale height along a circular sec- 
tion, and the cartesian coordinates (S, T, = z') are ob- 
tained from the system {x' , y' , z') by rotating around the 
z'-axis by an angle equal to (5* — 0), where ^' is the po- 
sition angle of the LM C bar major axis, roughly ~ 120° 
l(van der Marel II2OOII) . We fix it to 112.7°, so that the bar 
major axis S coincides with the minor axis of the disk ^. 

The column density, projected on the x — y plane is 
plotted in Fig. giving a global view of the LMC shape 
(disk and bar coplanar) for a terrestrial observer. We in- 
dicate the direction of the line of nodes, together with 
the positions of the microlensing events detected by the 
MACHO (filled stars and empty diamonds) and EROS 
(filled triangles)^ collaborations. The maximum value of 

^ The EROS group has publish ed the detection of 5 mi- 
crolensing events towards the LMC iLasserre et al. l200Ci) . The 
event EROS2-LMC-5 is dubious as its light curve is asymmet- 
ric. This candidate has actually been considered by MACHO 
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Fig. 1. Projection on the sky plane {x — y plane) of the 
column density of the LMC disk and bar. The numer- 
ical values on the contours are in lO^MQkpc"^ units. 
The three innermost contours correspond to 10, 20 and 
30 X lO'^M0kpc"^ The locations of the MACHO (black 
stars and empty diamonds) and EROS (triangles) mi- 
crolensing candidates are also shown. 

the column density, 41.5 x lO^MQkpc^^, is assumed in 
the center of the LMC. The positions of the MACHO mi- 
crolensing events in this x—y reference system are reported 
in Tabled 

In order to explore the consequences of different bar 
geometries, besides fro m the coplanar configuration a s in 
iGvuk et al I l|20nnl) and lWeinberg fc NikolaeTI l|200lh . we 
have considered two other possible geometries. In particu- 
lar, we drop the hypothesis that the disk and bar compo- 
nents are dynamically connect ed. In a first ca s e, dra wing 
inspiration from the paper of IZhao fc Evans"! l)2000() . we 
rotate the bar around an axis through the center, orthog- 
onal to the plane defined by the bar axis and the line of 
sight. As a second case we shift t he bar towards t he ob - 
server along the line of sight as in iNikolaev et al. 1 1)2004|) . 
In both cases we change accordingly the scale parameters 
of the bar so as to keep its observed projected shape on 
the sky plane fixed. For both configurations we considered 
several values for the shift and the rotation angle in or- 
der to study their influence on the microlensing quantities. 
In the following, as an illustration, we will give the results 
for two somewhat extreme cases, so that we will get upper 
limits for the corresponding microlensing quantities 

— a translation towards the observer of ~ 0.5 kpc; 



as it s candidate LMC-26 and rejected as being a supernova 
iMTlsztain .120031 . 



— a rotation of 40° with the south-east side towards the 
observer. In this case the linearly growing distance be- 
tween bar and disk is ^ 1.0 kpc at one scale length 
from the center. 

For the velocity dispersion of the stellar population 
of the bar, for which very few data are available, we 
take into account the qualitative results from numerical 
simulation that show a general trend of a higher line 
of sig ht velocity dispersion in th e central region of the 
LMC l|Zasov fc Khoperskov 12001 . We consider again two 
somewhat extreme cases: a line of sight velocity dispersion 
equal to that of the disk, dbjos = 20.2 km/s, and a second 
case with a higher value, CTb.ios = 30 km/s. 

3.3. LMC Halo 

We use two different models to describe the halo profile 
density: a spherical halo (model S) and an ellipsoidal halo 
(model E). The values of the parameters have been chosen 
so that the models have roughly the same mass within the 
same radius. 

• Model S. In this model we neglect the tidal effects 
due to our Galaxy, and we adopt a classical pseudo- 
isothermal spherical density profile, which allows us to 
compare the results obtained in this paper with the 
previous ones: 



Ph,s = Po,s 



1 



i?2 

a2 



OiRt-R), 



(3) 



where a is the LMC halo core radius, po.s the central 
density, Rt a cutoff radius an d t he Heaviside step 
function. We use a = 2 kpc (fAlcock et al. Il2000a() . We 
fix the value for the mass of the halo within a radius 
of 8.9 kpc equal to 5.5 x IO^Mq ijvan dcr Marcl ct al. 
I2OO2I) that implies /9o,s equal to 1.76 x 10^ M© kpc"^. 
Assuming a halo t runcation radius, Rt — 15 kpc 
ifvan dcr Marel et al. 2002), the total mass of the halo 
is « 1.08 X IQI^Mq. 

Model E. In this model the MW tidal effects induce a 
oblate spheroidal shape to the halo. We describe the 
halo density profile as 



Ph,E — Po,E 



1 



1 



9? 



(4) 



with a = 2 kpc, and the truncation radius Rt = 15 
kpc. For a LMC dark mass equal to 5.5 x 10^ Mq 
within a radius of 8.9 kpc, the central density is 2.20 x 
10^ M0 kpc"^, so that the LMC dark mass within the 
truncation radius turns out to be « 1.06 x 10^° Mq. 



3.4. Galactic halo 



We will consider for the galactic halo a spherical model 
with density profile given by: 



Pgh = Po 



Rl 



Rl 



(5) 
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where i? is the distance from the galactic center, Rc = 
5.6 kpc l|Primack et al. ]|T988) is the core radius, Rg — 
8.5 kpc is the distance of the sun from the galactic center 
and po = 7.9 x 10^ Mq kpc"^ is the mass density in the 
solar neighbourhood. 



4. Calculation of microlensing parameters 

As outlined in the introduction, here and in the following 
sections we want to fully exploit all the consequences of 
the geometry of the LMC, with the aim to shed new light 
on the still puzzling nature and location of the lenses de- 
tected in the microlensing surveys. In particular we study 
here two microlensing quantities, the optical depth and 
the microlensing rate. 

In our analysis, whenever we need to compare models 
and predictions with observational results, we arc going 
to use those presented by the MACHO collaboration only. 
The main reason is that this team has provided a complete 
description of their microlensing detection efficiency. 

From the total set of 17 MACHO microlensing events'^ 
one must exclude the event LMC-22, that is very likely a 
supernova of long duration or a n active galactic nucl eus 
in a galaxy at red-shift z = 0.23 l|Alcock et al. lEoOlcj) . 

We have detailed information on LMC-5: it is due to 
a lens located in the Galactic disk. Indeed t he lens proper 
motion has been observed with the HST ijAlcock et al. I 
l2001a|^ and the lens mass determined to be either ~ 0.04 
Mq or in the range 0.095 - 0.13 Mq, so that it is a true 
brown dwarf or a M4-5V spectral type low mass star. 
The other stars which have been microlensed were also 
observed by the HST, but no other lenses have been de- 
tec ted. This result is a lso co nfirmed by the analysis made 
by Ivon Hippel et al. 1 l)2003j) , by using optical and near 
infrared photometry on a subset of five lensing sources 
which are LMC main sequence stars or slightly evolved 
subgiants, (LMC-4, LMC-6, LMC-8, LMC-9, LMC-14). 
Their analysis rules out for main sequence lens masses 
> 0.1 Mq for distances out to 4 kpc. 

The events LMC-9 and LMC-14 are known to be due 
to lenses belonging to the LMC itself, i.e. to the bar, disk 
or halo component. The latt er event has been reco gnized 
thanks to the double source llAlcock et al. Il2001bl) . while 
the first shows the characteristic caustic crossing signature 
of a double lens I Alcock et al. lEoOObl) . For this reason 
we exclude the event LMC-9 from the following analysis, 
because we are interested in the study of an homogeneous 
set of single lens events. 

We remark that LMC-5, LMC-9 and LMC-14 are the 
only events for which it has been possible to make a deter- 
mination of the location of the lens. In Fig. ^ the events 
LMC-5 and LMC-9 are indicated with a special symbol, 
an empty diamond, whereas LMC-22 is not reported in 
the plot. 



^ The MACHO team actually claims the detection of 13-17 
events according to the chosen selection criteria. In the follow- 
ing we shall consider for our purposes the largest set. 
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Fig. 2. Plot of the interpolating function £{i) of the 
MACHO detection efficiency. 

We use the interpolating function £{t) of the mi- 
crolensing detection efficiency, calculated by the MACHO 
collaboration, as a function of event timescale t in th e in- 
terval 1 - 1000 days (Fig. 5 of I Alcock etall l)2000a|) ). It 
is shown in Fig. |21 by the continuous line, together with 
some points of the MACHO efficiency, for comparison. Let 
us recall that the MACHO definition of the duration time 
t is twice the Einstein time Te, the parameter we use in 
this paper. We get the following expression for £{t): 



£{i) ^ h{i) 6i[2.3 - Logt] ^[Logf] + 

h{i) 6'[Log£- 2.3] 6l[2.85 - Logi] 

/3(^>[Log^-2.85]^^[3-Log^] , 

with 9 the Heaviside function, and 



f2{i) 



Log' i 

£1.97 



iLog^t + ^Log^t 



/i(10^-3) 



(Logi - 2.3) 



-|^Q(2.3-Logt) 

h{i) = 0.06 + [/2(102 «5) _ 0.06] eio(2.85-Logt) 



(6) 

(7) 

(8) 

(9) 



As described in the previous section, the LMC has 
its own halo, with a tidal radius i?t = 15.0 ± 4.5 kpc. 
Assuming that the galactic halo extends to at least 65 
kpc, a fraction of the LMC halo total mass, enclosed 
in a sphere of radius Rt, is attributable to the galactic 
halo. A simple calculation gives the galactic halo mass en- 
closed in a sphere centered on the LMC with a radius 
of 8.9 kpc to be 9.7 x 10* Mq, a value corresponding 
to a sensible fraction of the LMC halo mass, of the or- 
der of ~ 17%. We will take into account this fact, by 
properly correcting the value of the central density of 
the LMC halo: for the spherical model the value is de- 
creased to po,s = 1-45 X 10''MQkpc~^ and the corre- 
sponding total LMC halo mass inside the tidal radius Rt 
is ~ 8.86 X 10^ Mq; for the ellipsoidal model the value is 
decreased to po,E = 1-81 x 10''MQkpc^'^ and the corre- 
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spending total LMC halo mass inside the tidal radius i?t 
is « 8.7 X lO^M©. 

4.1. Optical depth 

We now discuss the results obtained for the optical depth. 
The computation is made by weighting the optical depth 
with respect to the distribution of the source stars along 
the line of view. We use Eq. (7), Sect. 3.1, of Paper I, that 
we report here for self consistency reasons: 



3J6 -2.4 -12 1.2 2.4 3.( 



AttG 



poo rl 

Jo Jo 



D^l{Da,-Dal) 
Das 



Pi dDol 



■,dD„ 



(10) 



where pi denotes the mass density of the lenses, ps the 
mass density of the sources, D^i and Dos, respectively, the 
distance observer-lens and observer-source. 

The contour maps reported in Figs. 01 E [3 and refer 
to the case of coplanarity between bar and disk. In Fig. 
131 is reported the optical depth contour map for lenses in 
the Galactic halo, in the hypothesis that all the Galactic 
dark halo consists of compact lenses, together with the 
positions of the MACHO fields (see Fig. [T3| for the nu- 
meration of the fields), the microlensing events and the 
van der Marel line of nodes. We observe that almost all 
the fields (except three of them) fall between the contour 
lines corresponding to r = 46.3 x 10^^ and t = 49.3 x 10^*^. 
As expected the optical depth due to the galactic halo is 
a slowly variable function, and presents a slight near-far 
asymmetry: moving from the nearer to the farther fields 
along a line passing through the center and perpendicular 
to the line of nodes, the increase of the optical depth is of 
the order of « 6%. 

In Figs. 0] and El we report the optical depth contour 
maps for lenses belonging to the halo of the LMC, as- 
suming a spherical and an ellipsoidal shape, respectively, 
in the hypothesis that all the LMC dark halo consists of 
compact lenses. A striking feature of both maps is the 
strong near-far asymmetry. 

For the Model S, the maximum value of the optical 
depth, rniax,s — 8.05 x 10^^, is assumed in a point falling in 
the field number 13, belonging to the fourth quadrant, at 
a distance of ~ 1.27 kpc from the center. The value in the 
point symmetrical with respect to the center, belonging to 
the second quadrant and falling about at the upward left 
corner of the field 82, is Tg ~ 4.30 x 10~*. The increment 
of the optical depth is of the order of w 87%, moving from 
the nearer to the farther fields. 

The same happens for the model E: the maximum 
value of the optical depth, Tmax.E — 9.88 x 10~*, higher 
than the previous one, is found at about the same point, 
at the same distance from the center. In the symmetri- 
cal point with respect to the center, belonging to the field 
82, the optical depth is te — 5.05 x 10~^. The clUpsoidal 
shape of the LMC halo gives rise to a further enhancement 
of the near-far asymmetry, with an increase of the optical 
depth by « 95%. 




Fig. 3. Contour map of the optical depth for lenses in the 
galactic halo. The locations of the MACHO fields and of 
the microlensing candidates are also shown. The numerical 
values are in 10~* units. 



One can draw advantage from the different asymmetric 
behaviour of the optical depth in the two cases of lenses in 
the galactic halo or in the LMC halo, both to confirm the 
existence of a proper LMC halo and to disentangle the mi- 
crolensing events due to the Galactic halo from the ones 
due to the LMC halo. To this end, a good observation 
strategy would help, the goal being to allow the analy- 
sis of asymmetry of microlensing events belonging to two 
equivalent regions, placed symmetrically with respect to 
the line of nodes. An example of the kind of analysis and 
of results one could obtain is given in Sect.[Sl 

In Fig. El we report the optical depth contour map for 
self lensing, i. e. for events where both the sources and the 
lenses belong to the stellar population of the disk and/or 
the bar of the LMC. As expected, there is almost"* no 
near-far asymmetry and the maximum value of the optical 
depth, Tinax — 4.90 X 10^*, is reached in the center of 
the LMC. The optical depth then rapidly decreases, when 
moving, for instance, along a line going through the center 
and perpendicular to the minor axis of the elliptical disk, 
that coincides also with the major axis of the bar. In a 
range of about only 0.80 kpc the optical depth quickly falls 
to T ~ 2 X lO^'*, and afterwards it decreases slowly to lower 
values. 



* We note the peculiar "ear" shape of the contour line for 
T = 2 X 10~* due to the disk flaring. 
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Fig. 4. Model S: contour map of the optical depth for 
lenses in the LMC halo. The locations of the MACHO 
fields and of the microlensing candidates are also shown. 
The numerical values are in 10~* units. 

The calculated value of the optical depth in the center 
seems at first glance to be in contradiction with the value 
one gets using the Gould inequality ( Could ..199511) : 

r < 2 sec^ i ^ 1.34 x 10"^. (11) 

In fact, this inequality is derived making some simplifying 
assumptions which do not fully apply to our choice of the 
density profiles for the disk and the bar. In Sect.|H|we shall 
derive an improved version of the Could inequality, which 
will also apply to our density profile and leads to r < 
6 X 10^^, which does not contradict our above estimated 
value. 

Furthermore, we recall that the Could inequality is 
obtained under the hypothesis that the LMC is a virializcd 
system, that quite likely is note the case for some of its 
components. This should of course be taken into account 
when using it in comparison with the observations. 

As expected, the changes of the geometry as discussed 
in Sect. where we allow for a non coplanar morphol- 
ogy of the bar with respect to the disk, enhance the self- 
lensing optical depth in the bar region considerably (up 
to ~ 50%). On the other hand, the changes for the optical 
depth for lenses in the MW halo are negligible (at maxi- 
mum ^ 1%). For lenses in the LMC halo the variations, 
in the innermost region of the LMC, can be rather large 
(up to ~ 20%). However, the main feature that interests 
us, the near-far asymmetry due to the disk inclination, is 
not altered. Besides, as we are here mainly concerned with 
the self-lensing case, we do not discuss them any further. 



-3.6 -2.4 -1.2 12 2.4 3.6 




-3.6 -2.4 -1.2 12 2.4 3.6 



Fig. 5. Model E: contour map of the optical depth for 
lenses in the LMC halo. The locations of the MACHO 
fields and of the microlensing candidates are also shown. 
The numerical values are in 10~^ units. 



-3J6 -2.4 -1.2 1.2 2.4 3.6 




-3J6 -2.4 -1.2 1.2 2.4 3.6 



Fig. 6. Contour map of the optical depth for self Icnsing. 
The locations of the MACHO fields and of the microlens- 
ing candidates are also shown. The numerical values are 
in 10~^ units. The innermost contours correspond to val- 
ues 2.4 X 10-8, 3.2 X IQ-^, 4.0 x 10"^ and 4.6 x 10"* 
respectively. 
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Fig. 7. Contour map of the difference between the self- 
lensing optical depth calculated in the case of +40° ro- 
tated bar and the one with coplanar bar and disk. Going 
from the outer to the inner part, the contours correspond 
to increasing values of the difference: 0.1, 0.6, l.landl.6x 
10~*. The difference is null outside the bar region. 



Fig. 8. Contour map of the difference between the self- 
lensing optical depth calculated in the case of 0.5 kpc 
shifted bar and the one with coplanar bar and disk. Going 
from the outer to the inner part, the black contours cor- 
respond to increasing positive values of the difference: 
+0.1, +0.5, +0.9, and +1.3 X 10"^; the lighter contours to 
negative values of the difference: —0.1, and — 0.2 x 10^^. 



In any case, outside the bar region, these differences drop 
abruptly to zero. 

For a rotated bar. Fig. [7| the increment of the self- 
lensing optical depth, with respect to the symmetric copla- 
nar case, is slightly asymmetric with respect to the bar 
major axis. This is due to the different variations of the 
source-lens distances between the west side, where the 
sources are in the bar and the lenses in the disk, and the 
east side, where the opposite happens. The relative incre- 
ment can be as large as ~ 50 %. 

When we allow for a translation of the bar, Fig.|Sl the 
variation of the optical depth, with respect to the copla- 
nar case, is higher in the region below the bar major axis 
(south-west side), even if the projected density of the bar 
is the same as shown in Fig. ^ This is induced by the 
inclination of the disk plane, that gives rise to increasing 
distances between the bar and the disk plane for lines of 
sight towards the lower part of the bar. In this case we find 
also a small region with a negative difference of the opti- 
cal depth, shown by the lighter contour lines (north-east 
side). The relative increment can be as large as ~ 35%. 

It is interesting to note that a larger statistics of ob- 
served events in the central region might eventually allow 
one to discriminate between the different bar models. 

We notice also that, due to the loss of the symmetry, 
the Gould inequality no longer applies in these cases. 



4.2. Self tensing microlensing event rate 

The distribution , the differential rate of microlensing 
events with respect to the Einstein time T^, allows one to 
estimate the expected typical duration and the expected 
number of the microlensing events. In this paper we evalu- 
ate the microlensing rate in the self-lensing configuration, 
i. e. lenses and sources both in the disk and/or in the 
bar of the LMC. We take into account also the transverse 
motion of the Sun and of the source stars. 

We assume that, to an observer comoving with the 
LMC center, the velocity distribution of the source stars 
and lenses have a Maxwellian profile^. Moreover, since 
the component parallel to the line of view is irrelevant to 
microlensing, we integrate over the parallel component. 
Therefore, the two-dimensional transverse velocity distri- 
bution of the source and lens stars is: 

1 (-^) 

/(t3_L) d^^^L = ^'lasJ v^dv^dip , (12) 

^ The Maxwellian profile of the velocity distribution is the 
first term of a series expansion in terms of Gauss-He rmite 
moments (Ivan der Marel fc Frannil993l: iGerhard Ill993ll . See 
Section 3.2 of Paper I. 
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where we have used cylindrical coordinates^ d'^v± — 
v± dv± d(j). 

With respect to an observer comoving with the LMC 
center, the transverse velocity of the microlensing tube 
at position Z?oi ~ x Dos (0 < x < 1) is given by 
Vt± = {1 — x)vqj_+xVs±, where is the velocity of the 
source, and vqj_ represents the Sun velocity in the plane 
orthogonal to the line of sight, as measured by an observer 
comov ing with the LMC center. It results I'WqxI — 406 
km/s l|van der Marel et al. llioO^ . 

Let us consider a segment of cylindrical ring at position 
x, of length dDoi, with radius equal to the Einstein radius 



Re 







flDosX{l 



angular amplitude da and thickness {vii_ ■ n) dt = 
vi± cos 6 dt, where — f < ^ < f is the angle between 
the inner normal to the cylindrical ring segment and the 
vector vij_. v\j_ represents the lens velocity in the plane or- 
thogonal to the line of sight, as measured by an observer 
comoving with the microlensing tube at position x. We 
use solar mass units, defining fi — jj"; where M is the 

lens mass. The number of lenses with mass in the range 
(/X, n + d/i) inside the volume element of the cylindrical 
ring segment is equal to the product of the differential 
number density of the lenses d /i by the volume ele- 
ment vi± cos Re dt da dDoi- The fraction entering in the 
time interval dt in the microlensing tube with transverse 
velocity vi±, is given by those lenses that with respect to 
the LMC reference frame have velocity 



*1_L = Vi±_ + Vtl_ = Vll_ + x) Vqj^ + X Vsl_ 



(13) 



In the self-lensing case, a; ~ 1, and to first approxima- 
tion we can neglect the term proportional to (1 — x). As a 
consequence the velocity distribution with respect to the 
microlensing tube, in the general case of two different val- 
ues for the velocity dispersion, is given by: 



vij_ dv\± d0 



2 IT [pd{x) + Phix)] 
Pdjx) 

vi± dvi± dO 
2tt [pd{x) + ph{x) 

where ad = crd,ios: 



Ph{x) 



■h{pd,crd] Pb,o-b; x,vij_,Vs±; a, 6) (14) 



CTb = cTbjos and a is the angle that 
the segment of cylindrical ring at position x, of angular 
amplitude da, forms with the direction of the orthogonal 
component to the line of sight of the velocity vector of the 
source star, Vs±. 

We need now to specify the form of the number den- 
sity. Assuming that the mass distribution of the lenses is 
independent of their position in the LMC {factorization 

^ Eq. 1121 gives the correct mean square transverse velocity 
< vl >= 2crL. 



hypothesis l|de Ruiula et ainil995|) ). the lens number den- 
sity per unit mass is given by 



dn\ 
dp 



Pd 



M 







Ph dnp 
dp 



(15) 



where we use as given in IChabrier I l)200l[l . We con- 
sider both the power law and the exponential initial mass 
functions^. However, we find that our results do not de- 
pend strongly on that choice and hereafter, we will discuss 
the results we obtain by using the exponential IMF only. 

Let us note that, considering the experimental condi- 
tions for the observations of the MACHO team, we use as 
range of variability for the lens masses 0.08 < p < 1.5. 
Namely the lower limit is fixed by the requirement that 
the lenses belong to the luminous component of the LMC, 
while the upper limit is fixed by the requirement that the 
lenses are not resolved stars^. 

The total differential rate dF at which lenses enter 
the portion of the microlensing tu be at position along 
a fixed fine of sight, is given by ijde Riiiula et al. Ill99lt 
iGriest II199I : 



dF 



vf± dvij_ 
47rWM0 Jo 



dip / d9 cost 



Dos dDos I dx I dp -J— Re 



1.5 

i„: 

0.08 d/i 



dwsj_{i. 



PdjDos) -f^ , PhjDos) -t^ 

^ — e d -I T. — e b 



^2 7r 

/ daft.(pd,crd; Pb,crb; a;,ui_L,Us_L; a, 
Jo 



0) (16) 



where 



Af^ / dDos [pdiDos) + PhiDos)] , 

having assumed that the number of detectable stars varies 
with the distance as D~^. The integration limits dmin, 
rfmax represent the distances from the observer of the in- 
tersection points of the line of sight with the LMC tidal 
surface. 

Finally, as we are interested in the distribution 
we change variable from vi± to Tg, bearing in mind that 
vi±^ — After integration on ip, 9 and a, and taking 
into account the detection efficiency function, Eq. ©, we 
obtain: 



dF 
dJ^ 



dF 32G^Mq£{Te) 
dTE 



dDosDts 



^ We have used the same normalization as in Paper I with 
the mass varying in the range 0.08 to 10 Mq. 

^ We have checked that the results are insensitive to the pre- 
cise upper limit value. This is also confirmed by the following 
analysis on the expected value for the mass of the lenses. 
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1 rl.5 
2 /-I „a2 



L 



dxx (1 — x) I dfifi 

0.08 "^/^ 



2 dno 



Pb(a;)pb(£'os) 
— ^ -t: — TT- e ^ 



2 G Mq Dos a; (1-a;) fi 



/ 2 G 

-T^ 7^ e V d 



2 G Dos X (1-a;) 



^2 + x^al 



Ph 



/ 2 G 



■2+x2„2) ^2 



2 G Dos a: (1-x) ^ 



Pd(a;)Pd(-Dos) 



2 G Mq Dob 3: (I-k) 



(17) 



In the 
tions 

/ d^r 



(i + x)V2 

following section we need also the two distribu- 



\ dTE da; 

/ d^r 

\ dTs d^ 



d^r 



dTs dx 



d^r 



dTEd// 



(18) 



(19) 



calculated assigning to Te the corresponding effective 
measured value. 



4.3. A statistical analysis for self-lensing events 
discrimination 

In this section we will show that, in the framework of 
the LMC geometrical structure and dynamics outlined in 
Sect. 13 a suitable statistical analysis allows us to exclude 
from the self-lensing population a large subset of the de- 
tected events. To this purpose, assuming all the 15 events 
as self lensing, we study the scatter plots correlating the 
self lensing expected values of some meaningful microlens- 
ing variables with the measured Einstein time or with the 
self-lensing optical depth. The idea underlying this analy- 
sis is based on the search of self-lensing average collective 
properties with different behaviour in the two different re- 
gions of the LMC: the bar with its nearby neighbourhood 
and the disk region external to it. In this way we even- 
tually show that a large subset of events is incompatible 
with the self-lensing hypothesis. 

We have calculated the self-lensing distributions 
of the rate of microlensing events with respect to 

the Einstein time Tg, along the lines of sight towards the 
15 events of the MACHO collaboration, in the case of a 
Chabrier exponential type IMF. As an example we show in 

Fig- El the distributions (^g^^ calculated along the lines 

of sight pointing towards the events LMC-1 (solid line), 
LMC-8 (dashed hue) and LMC-23 (dot dashed hne)^. 

® We have also used the Chabrier power law IMF, obtaining 
slightly higher values, of the order of 10%, with respect to the 
exponential one. In the following we have used constantly the 
Chabrier exponential IMF. 




50 100 150 200 S50 300 350 400 

Te Cdays) 

Fig. 9. Differential rate of the microlensing events with 
respect to the Einstein time Te, along the lines of sight 
pointing towards the events LMC-8 (dashed red line), 
LMC-1 (solid black line) and LMC-23 (dot dashed blue 
line). The y-axis values are in 10"^^ units. 



With these distributions we have calculated the me- 
dian^° Te^ so % and the values Te^ ig % and Te^ §4 % that sin- 
gle out the extremes of the 68% probability range around 
the median (not to be confused with a la error). In Table 
^we report these values for each observed MACHO event. 

Trying to see if the geometry can help, we have ana- 
lyzed how the self-lensing expected values of Te depend 
on the position, or better still on the optical depth, taking 
into account that the LMC disk symmetry is elliptical and 
not circular. 

In Fig. ^] we report on the y-axis the observed values 
of Te (empty boxes) as well as the expected values for 
self lensing of the median Te^5o% (filled circle) evaluated 
along the directions of the events. On the x-axis we re- 
port the value of the self-lensing optical depth calculated 
towards the event position. The optical depth increases as 
one moves from the outer regions towards the center of 
the LMC according to the contour lines shown in Fig. El 
An interesting feature emerging clearly is the decreasing 
trend of the expected values of the median Te^5o%, going 
from the outside fields with low values of tsl towards the 
central fields with higher values of rgL. The variation of 
the stellar number density and the flaring of the LMC disk 
certainly contributes to explain this behaviour. 

We now tentatively identify two subsets of events: the 
ten falling outside the contour line tsl = 2 x 10~^ of Fig. 
[Bland the five falling inside. In the framework of van der 
Marel et al. LMC geometry, this contour line includes al- 
most fully the LMC bar and two ear shaped inner regions 
of the disk, where we expect self-lensing events to be lo- 
cated with higher probability. 

At first glance, we note that the two clusters have 
a clear-cut different collective behaviour: the measured 



As the distribution is strongly asymmetric, to describe the 
expected value of Te we use the median value, a more mean- 
ingful parameter than the average value. 
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event 


X (kpc) 


y (kpc) 


TE.obs 


2e,50% 


rE,i6% 


2e,84% 


1 


1.017 


0.909 


22.3 


64 


33 


126 


4 


0.746 


-0.814 


29.5 


66 


35 


128 


5 


0.797 


-0.559 


49.1 


65 


34 


127 


6 


0.102 


-0.423 


59.5 


55 


28 


111 


7 


1.796 


0.189 


66.8 


73 


38 


140 


8 


0.185 


0.062 


43.1 


48 


25 


95 


13 


0.280 


0.916 


65.0 


66 


33 


128 


14 


-0.523 


-0.487 


65.0 


51 


26 


103 


15 


1.652 


0.048 


23.9 


72 


38 


137 


18 


-1.253 


-1.168 


48.2 


72 


38 


138 


20 


2.478 


-0.316 


47.2 


77 


40 


146 


21 


2.324 


0.211 


60.5 


77 


40 


146 


23 


1.517 


-1.037 


55.4 


72 


38 


138 


25 


2.072 


1.517 


55.4 


79 


41 


149 


27 


1.619 


0.388 


32.8 


67 


34 


130 



Table 1. Microlensing rate: the results for the predicted Einstein time. Case with disk and bar coplanar and bar 
velocity dispersion <t\, — 20.2 km/s. For each observed MACHO event we report for the Einstein time (days): the 
observed value and the predicted median value with the two extreme of a 68% probability range around it. In the 
second and third column we report the position in the reference frame centered in the LMC as defined in Sect.|21 



Einstein times of the first 10 points fluctuate around a 
median value of 47 days, very far from the expected values 
of the median Te, ranging from 65 days to 79 days, with 
an average value of 72 days. On the contrary, for the last 
5 points, the measured Einstein times fluctuate around a 
median value of 51 days, near to the average value 57 days 
of the expected medians, ranging from 48 days to 66 days. 
Let us note, also, the somewhat peculiar position of the 
event LMC-1, with a very low value of the observed Te. In 
the following analysis it will be shown that most probably 
this event is homogeneous to the set at left of the vertical 
line in Fig. 1101 and it has to be included in that cluster. 

This plot gives a first clear evidence that, in the frame- 
work of van der Marel et al. LMC geometry, the self- 
lensing events have to be searched among the cluster of 
events with tsl > 2 x 10^®, and at the same time that 
the cluster of the 10 events plus LMC-1 belongs, very 
probably, to a different population. 

In order to further improve our statistical analysis, 
we have calculated the distributions (^g^^^^ along the 

lines of sight pointing towards the 14 LMC microlens- 
ing events^^, taking, for each line of sight, the observed 
Einstein time value. As an example we show in Fig. ^] 
this distribution calculated for the two events LMC-6 and 
LMC-4. LMC-6 is representative of the events for which 
the modal value /imod falls inside the self-lensing mass in- 
terval 0.08-1.5. The second one has been chosen to put 
in evidence that there are also events for which /imod falls 
outside, in the range 0-0.08. 

Fig. El is the scatter plot between the measured Te 
(x-axis) and the modal value of the lens mass /i = 

(y-axis) , calculated by the distribution ( ^ji^^^ ) . In this 




The previous 15 events minus the event LMC-5, formerly 
recognized as a Galactic disk event. 



Fig. 10. Scatter plot of the observed (empty boxes) val- 
ues of the Einstein time and of the expected values of 
the median Te^5o% (filled stars), with respect to the self- 
lensing optical depth evaluated along the directions of the 
events. The Te^5o% for the galactic disk event LMC-5 is 
represented by a filled box. 

case we prefer to use the modal value rather than the 
median one, because it is independent of the amplitude 
of the interval of the allowed values of the lenses mass, 
whereas the median value varies in accordance with this 
choice. 

We find that the events LMC-1, LMC-4, LMC-15 and 
LMC-27 have a modal value of the lens mass smaller than 
the lower limit. We consider this result as a strong indi- 
cation to exclude these four events from the self-lensing 
class. We have then calculated the linear correlation be- 
tween ^rnod and Te for the 6 remaining points out of the 
cluster of ten. We find a high linear correlation, as shown 
by the dashed green straight line in Fig. El and by the 
calculated linear correlation coefficient equal to 0.963. We 



12 



L. Mancini et al.: LMC Self-lensing from a new perspective 




Fig. 11. Differential rate of the microlensing events with 
respect to the mass of the lens in direction of events LMC- 
4 and LMC-6. The dashed vertical line is at ^ = 0.08. The 
y axis values are in 10^'^ units. 
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Fig. 12. Scatter plot of the modal value of fi (y-axis) with 
respect to the measured Einstein time (x-axis) of each 
event. The label identifies the event. Filled (empty) trian- 
gles represent points with tsl > 2 x 10^^ (< 2 x 10^^). 
The horizontal line indicates the lower limit of the lens 
mass for self lensing, fi = 0.08. The four events LMC-1, 
LMC-4, LMC-15 and LMC— 27 have a modal value of 
the lens mass smaller than the lower limit and therefore 
are not represented. The dashed line represents the corre- 
lation line for the 6 events with tsl < 2 x 10~*. 



observe that the values of fJ-mod of the six events range be- 
tween 0.08 and 0.13, an interval narrow enough to justify 
the linear approximation to represent a small portion of 
a parabolic curve, around which we expect the correlated 
points to disperse, bearing in mind that Te is proportional 
to y/JI. The six linearly correlated events, therefore, con- 
stitute a homogeneous population, clearly distinct, also in 
this parameter space, from that formed by events with 
TSL > 2 X 10~^ for which, for a given value of the ob- 
served Einstein time, we get significantly higher values for 
the mass. 

We have calculated also the distributions ^ ^j-Jdx ) 
along the lines of sight pointing towards the 14 LMC mi- 



0.03 
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0.015 



0.01 




10 50 

TE.pte fdays) 

Fig. 13. Correlation plot of the median value (square 
boxes) oi 6 — (y-axis) with respect to the mea- 

sured Einstein time (x-axis) of each event. The label iden- 
tifies the event. The dashed line is the correlation straight 
line of the cluster of 10 events. 



crolensing events, in order to obtain a second independent 
check of the homogeneity of the cluster of events with 
< 2 X 10~®. Fig. ^1 is the scatter plot between the 
measured Te (x-axis) and the median (square boxes) of 
the parameter S = — l — ^ (y-axis), proportional 

to the distance between the lens and the source. 

We find that the points of the first cluster have a high 
linear correlation, as shown by the calculated linear cor- 
relation coefficient equal to 0.955. We observe also that 
adding the event LMC-1 to the set of 9 and recalculat- 
ing the linear correlation coefficient between S and T^, 
we find a small increase, from 0.955 to 0.965, suggesting 
that this event, lying inside the contour line at optical 
depth TgL = 2 X 10^*, is very probably homogeneous to 
the population of the 9 falling in the outer part. The ten 
events, represented by gray square boxes, are highly cor- 
related as shown by the dashed green straight line in Fig. 
1131 This is, again, a strong indication that they constitute 
an homogeneous population of events. Together with the 
fact that, as shown in Fig. 1101 the measured Einstein time 
fiuctuates around a median value very far from the median 
Einstein times, calculated with the self-lensing formulae, 
this allows us to exclude that this class of events belong to 
self lensing. But, before any definitive assessment on the 
nature of these events can be made, we have to wait for 
an analogous statistical analysis for the case of microlens- 
ing events due to lenses in the halos of the LMC or of 
our Galaxy: such an analysis is now in preparation for a 
forthcoming paper. 

We come now to the discussion of the changes in the 
microlensing rate induced by the different bar geometry 
configurations and velocity dispersions, as introduced in 
Sect. 13.21 The main point to be stressed is that the sepa- 
ration in the two populations for the events already iden- 
tified is in these cases enhanced. Indeed, the expected 
characteristics change significantly only for those events 
along the lines of sight pointing towards the central region: 
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event 


TE.obs 


Te,50% 


2e,50% 


2e,50% 


r(10-») 


r(10-») 


r(10-») 






coplanar 


rotated 


shifted 


coplanar 


rotated 


shifted 


1 


22.3 


64 56 


76 67 


64 56 


2.24 


3.33 


2.23 


6 


59.5 


55 47 


57 48 


62 53 


2.84 


3.16 


3.90 


8 


43.1 


48 37 


53 41 


51 40 


4.72 


5.99 


5.58 


14 


65.0 


51 41 


60 47 


56 44 


3.72 


5.33 


4.53 


27 


32.8 


67 59 


75 67 


69 61 


1.75 


2.29 


1.94 



Table 2. For each MACHO event located within or nearby the bar region, vi^here the variation due to the change in the 
geometry turn out to be significant, we report for the Einstein time (days): the observed value and the predicted values 
for the three geometry configurations, bar coplanar with the disk, bar rotated by 40° with respect to the disk plane 
and bar shifted by 0.5 kpc along the line of sight towards the observer. The left values are for a velocity dispersion 
of the bar component of (Tb=20.2 km/s (as for the disk), whereas the right values are for = 30.0 km/s. We report 
also the values for the self-lensing optical depth in the direction of the events. 



LMC~6, LMC-8 and LMC-14, and, marginally, LMC-1 
and LMC-27, as it is the case for the self-lensing optical 
depth. In Table El we report the predicted values for the 
Einstein time for these 5 events, the changes for the others 
being at most of ~ 5% (in particular this is the case for the 
event LMC-13). As expected, an increase in the velocity 
dispersion for the bar component leads to a decrease (^ 
20%) for the predicted values of the Einstein time. For the 
events LMC-6, LMC-8 and LMC-14 the predicted values 
just fluctuate around the observed ones, while they remain 
substantially different for the events LMC-1 and LMC-27. 
The change in geometry leads, on the contrary, to a rise 
for the expected values of the Einstein time, at most ~ 
10% (the effect being stronger for the rotated bar), for 
the events LMC-6, LMC-8 and LMC-14. For the event 
LMC-1 the difference between the observed and the pre- 
dicted value is further enhanced (up to ^ 16%). This gives 
further support to the guess that this short duration event 
is likely to belong to a different population. We note that 
the two opposite effects on the predicted Einstein time, 
the decrease linked to the rise of the velocity dispersion 
opposed to the increase linked to the non coplanar geom- 
etry, are actually to be expected as F cx r/T. Analogous 
considerations emerge from the analysis of the differen- 
tial rate with respect to x and the mass. In particular, 
for the mass, we note again two opposite effects. A rise 
induced by the increase in the velocity dispersion, and a 
decrease linked to the non coplanar geometry. We observe 
a maximum increase of ~ 15% for the event LMC-14. 

To complete this analysis let us make a further point. 
In the differential rate, Eq.EI we get an average result for 
all the possible configurations with lens and source either 
in the bar or in the disk of the LMC. This is of course 
coherent with our approach aimed at the study of the sta- 
tistical properties of the self lensing population as a all. 
However, it is in principle interesting to notice that the 
expected characteristics can be actually rather different 
depending on the configuration. As an example, we con- 
sider the event LMC-14, coplanar bar geometry, Ub — 20.2 
km/s. With respect to the previously obtained result, ex- 
pected median value for the Einstein time 7e^5o%=51 days, 
we get a shorter value in the lens-source bar-bar configu- 



ration, 43 days, a significantly longer value for the disk- 
disk configuration, ~ 67 days, and an intermediate one, 
55 days, for both mixed configurations, disk-bar and 
bar-disk. We recall that the observed value for LMC-14 is 
rE=65 days. 

4.4. Expected number of self-lensing events 

We compute the "fi eld exposure", E'ficid, defined, as in 
lAlcock et al. I l)2000aj) , as the product of the number of dis- 
tinct light curves per field by the relevant time span, pay- 
ing attention to eliminate the field overlaps. Furthermore 
we calculate the distribution along the line of sight 
pointing towards the center of each field. In this way we 
obtain the number of expected events for self lensing, field 
by field, given by 

iVsL,fiold = Afield / :J77^ f (Te) dTE , (20) 

Jo QTe 

where £{Te) is the detection efficiency. 

In Fig. 1141 we report in parenthesis the expected num- 
ber of events in each field we obtain by using the Chabrier 
exponential IMF. Summing over all fields we find that the 
expected total number of self-lensing events is ~ 1.2, while 
we would get ^1.3 with the the double power law IMF, in 
both cases altogether 1-2 events. Clearly, taking also into 
account the uncertainties in the parameter used following 
the van der Marel model for the LMC the actual num- 
ber could also be somewhat higher but hardly more than 
twice our estimate. A conclusion we had already reached 
in Paper I. 

We have also computed the influence of the bar geom- 
etry and dispersion velocity on the number of expected 
events. We find at most an increase of ~ 50%, bar rotated 
and (7h— 30 km/s, that gives altogether about 2-3 events. 

In conclusion, we have identified, according to their 
expected characteristics as compared to the observed one, 
a set of 4 possible events belonging to the self-lensing 
population (LMC-6, LMC-8, LMC-13 and LMC-14). As 
already observed, there are some experime ntal evidences 
that t he LMC-14 is a self-lensing event ("Alcoc k et al. I 
l2001bi) . If this is true, our prediction of at most 2-3 such 
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events seems thus to indicate that some of the others are 
more hkely not to be due to self lensing. 



5. Asymmetry 

In this section we consider from a different perspective the 
issue of the expected near/far asymmetry for events due to 
lenses located in the LMC or the MW halo. In particular 
we study, by means of a statistical analysis, the spatial 
distribution of the events, both in the framework of the 
van der Marel picture and, for comparison, of the older 
one. The aim here is to give further support to our main 
conclusion, namely that self lensing, for which we do not 
expect any asymmetry in the spatial distribution, can not 
explain all the detected events. 

In Fig. El the 30 well-sampled fields analyzed by the 
MACHO collaboration (red squares), together with the 16 
events^^, are located in a reference frame (black axes) cen- 
tered in acM = 5^27.6", Scm = - 69.87°, J2000. We also 
report the position (triangles) of the microlensing events 
found by the EROS collaboration (Lass_erre et al. 2Q00) 
even if we do not consider them in the present analysis. 

We divide the LMC field in three regions: a strip cen- 
tered on the line of nodes and delimited by two parallel 
straight lines at a distance of w 0.35 kpc from the line 
of nodes, and two outer regions, belonging to the closer 
north-east side and to the farther south-west one. The 
amplitude of the central strip refiects the uncertainties 
in the positio n of the center o f mass of carbon stars, as 
calculated bv lvan der Marefet al. . 1.2002l) . at one a level. 
The exclusion of a substantial part of the bar region im- 
plies that the following discussion is not affected by the 
different choices for the bar geometry. 

The green dashed-dot line represents the line of nodes 
with a position angle 8 = 129.9° and the two green dashed 
lines delimit the two outer regions belonging respectively 
to the near and far part of the LMC. 

On the same figure a second reference frame is also 
reported, according to pre-van der Marel LMC models. 
The frame (light blue a xes) is centered iii ofgas = 5^17.6" 
S = -69°02', J2000 (|Kim 



9981) . a nd the line of 
node s has a position angle Gcyuk = 170° 1 Cvuk et al. I 
region of exclusion, similar to the previous one 
will be considered also in this case, even if not drawn on 
the figure for clarity sake. 

At first glance one observes that the distribution of the 
events shows a clear near-far asymmetry in the post-van 
der Marel geometry, namely, they are concentrated along 
the extension of the bar and in the south-west side of 
the LMC. On the contrary, the asymmetry is almost com- 
pletely lost in the pre-van der Marel geometry, where the 
distribution reflects almost exactly the different weights of 
the observed fields in the two half planes, as we will show 
in the following. 



The little empty circle in the field number 5 locates the 
position of the baricenter of the 16 MACHO events, whose 
coordinates in the post-van der Marel reference frame are 
(0.96 kpc, - 0.02 kpc). 

5.1. Near- far asymmetry of the observed microlensing 
events in the LMC 

In order to give a quantitative analysis of the near-far 
asymmetry, we will concentrate on the MACHO observed 
fields (and the corresponding observed events) in the outer 
regions located on the opposite sides of the line of the 
nodes and external to the dashed lines, which, at a good 
confidence level, belong respectively to the near and far 
part of the galaxy. 

As a first point we determine the fraction of the 
MACHO fields included respectively in the north-east 
closer region "1" and in the south-west farther region "2" , 
and we calculate in each the quantities Ei and £2, defined 
as the sum of the product of star number pro field, times 
the corresponding observation time, where we count only 
once the stars in the intersecting part of the fields. The 
ratio I]i/(Si -I- S2) gives us the probability pi that a mi- 
crolensing event would fall in the first or second region, 
respectively. Let us note that the probability scheme so 
defined depends only on the global observation strategy 
in the near and far region (fields distribution and observa- 
tion time). We find pi « 0.31 in the case of the post- van 
der Marel LMC geometry, and pi « 0.48 for the pre-van 
der Marel LMC geometry. 

We are interested in testing whether the observed 
events support the modelled probability schemes. The 
Pearson chi-square statistic An provides a non parametric 
test for the comparison: 



{n, -Np,f 
Np^ 



(21) 



The ensemble of 17 events discussed in Sect. |1] with the 
exclusion of LMC-22. 



where rii are the events falling in the i-region, and N — 
ni + n2. An is approximately distributed as chi-square 
with one degree of freedom. Small values of An tend to 
support the null hypothesis that the pi match the mea- 
sured distribution. 

In the post-van der Marel geometry we get, respec- 
tively, 1 and 9 observed events in regions "1" and "2", 
and then 

Aio = 2.11. 

On the contrary, in the case of pre-van der Marel geometry, 
we obtain: 

A12 = 0.22, 

starting from a set of 5 and 7 observed events. 

In the pre-van der Marel geometry An gets a value 
near to zero. At the confidence level of ~ 64% the null 
hypothesis that the events distribution reflects almost ex- 
actly the weights of the two regions has to be accepted. 
This implies that the distribution of the lenses should be 
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Fig. 14. Locations of the 30 MACHO fieMs with respect to post- Van der Marel LMC geometry (black reference axes) 
and with respect to pre-Van der Marel LMC geometry (light blue reference axes) . Also shown are the locations of the 
MACHO and EROS microlensing candidates. The distances on the axes are in kpc. For each field, inside brackets, we 
report the expected number of self-lensing event as discussed in Sect. 14.41 



almost homogeneous along the lines of sight through the 
different regions of the LMC. 

On the contrary, in the post-van der Marel geometry 
A]^ assumes a value far enough from zero. At the confi- 
dence level of ~ 85% the null hypothesis must be rejected. 
This means that the observed asymmetry is greater than 
the one expected simply on the basis of the observational 
strategy. 



We are aware that these results have to be treated 
with caution inasmuch as the number of events is small. 
However, we note that the observed near-far asymmetry 
is coherent with that induced by the inclination of the 
LMC disk already discussed, looking at the optical depth, 
in Sect.Ol 
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6. An improved Gould inequality 

ICould I l(l995l) ingenious calculation provides two inequal- 
ities (Eq. (2.8) and Eq. (3.3) of his paper), the second 
one involving the Jeans equation and the virial theorem, 
which allow, in some cases, a quick evaluation of an up- 
per limit to the self-lensing optical depth along a line of 
view through the center of a galaxy. In particular t his was 
applie d to the LMC disk, in order to exclude the ISahu I 
l)l994() hypothesis that the observed optical depth could 
be fully explained by self lensing. 

Our aim here is to obtain an improved version of Gould 
inequality, and, at the same time, clarify its limits of ap- 
plicability. We start from Eq. 1)10(1 and consider the case 
of a line of view passing through the center of the LMC. 

Let us assume that the lens mass density p\ be a func- 
tion of an homogeneous polynomial of degree n\ in the vari- 
ables (^, 77, C) of the reference frame defined in Sect. 13.1 1 
and analogously that the star mass density ps be a func- 
tion of an homogeneous polynomial of degree Us^^ . Let us 
also denote with z the current coordinate along the line 
of view through the center of the LMC disk, which we 
assume as the origin of the z-coordinate, and with i the 
disk inclination angle. Keeping in mind that for points 
belonging to a line of view passing through the center it 
results: 



rj — z sin i 
C — z cos i 

we can write for the lens density distribution 



p;[0, -n, C] = Pi 



En 
k=0 



.(1) 



^k,n — k 



n 

')l,o 



Pi 



71 



with 

71 



k=0 ^^w^ 



sin'=(^) cos"-'^-(*) 



(22) 



(23) 



and rjifi, ^i^o the scale lengths of the lenses density distri- 
bution. Analogous expressions can be written for the star 
density distribution, changing properly the suffix "1" with 



From now on we follow iGouldl l)l995|) . and put Doi = 

d + z and Dos = d + w, where d is the distance, along the 
line of view, from the observer to the center of the disk. 
Taking into account that 

d + z 
d + w 

we obtain the following inequality for the optical depth: 

4-rG f* 



(24) 



/ dz (w — z) pi 

•J —d 








73 





This is implicit in the lGouldl ^99^ derivation. 



where t is the the LMC tidal radius. 

We now integrate by part twice, a first time obtaining: 

4ttG f* 



r<^/!,d^p. 



7 

71 ('0 



« (fi) 
(-)■" 



J_^dwp^ (^)' 

and a second time obtaining: 
AttG 



(25) 



T < 






7i(*) 




/ dzpi 




J-d 





I-d Pi 




ni- 


I-d Pi 


it;) 


rii- 


/r^dz 


Ps 




1 "a- 


I-t P- 




d 



(26) 



In case of self lensing ps and pi coincide. Moreover, for 
distances higher than the tidal radius of the LMC, the 
star mass density ps = 0, and we can move the lower 
integration limit from — < or —d to —00, and the upper 
from t to +00. We are considering a line of view passing 
through the center of the LMC, therefore we can use the 
symmetry property of the mass density distributions with 
respect to change of sign of the parameter z and we obtain: 



— J— 00 







7s(*) 


(&)1 



where 



Fs{v) 



2j^dzps 




lis- 









(27) 



(28) 



A suitable change of variables gives: 
2 7rGC.^n 



T < 



where 



Fs{v) 



[ls{^)]-- ■ 
Ps(u"=) du / [l~F^{v)\ dv, 



2/;p, (u"-)du 
!Too ps du 



(29) 



(30) 



and the integrations are now made on adimensional vari- 
ables. 

Let us observe that the inequality (|29|) can be applied 
for any inclinat ion angl e i of the disk plane, not only for 
i « 7r/2, as in iCould I ljl994 - namely, we have no diver- 
gence problem for i tt/2. For instance, in the case of 
a double exponential disk, with scale lengths respectively 
equal to i?d and Cd, we obtain: 



27rG<:i 
5 



[cos(j) + -^sin(i)] 

J^-[l-FHv)] dv. 



(31) 
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Or in the case of a gaussian profi le of the mass densi ty dis- 
tribution for the LMC bar, as in ICxvuk et al. I ll200(t . with 
scale length rj^, along the axis bar and Cb in the orthogonal 
section, we obtain: 



sin^(i) 



(32) 



For a bar h aving boxy contours a nd sharp edges as in 
the paper of IZhao fc Evansi l)2000() . parameterized by an 
exponential profile with a fourth degree polynomial at the 
exponent, with scale length 77b along the axis bar and 
in the orthogonal section, we obtain: 



^ T 



cos''(i)+(^)%in''(i) 



d^^- (33) 

Eg. 129|l is the improved version of Eq. (2.8) of lGouldl 
lll9951. 



Proceeding in the same way as in iGoiild I l)l995|) . we 
obtain the second Gould inequality, relating the optical 
depth with the mass-weighted velocity dispersion. 



r<2^-^[7s(z)]-- 



(34) 



Let us note that these inequalities are based on two 
properties: 

— that the star density be a function of an homogeneous 
polynomial of the variables, of any degree; 

— that the line of view pass through the center, in such 
a way we can use the symmetry by reflection of the 
mass density distributions. 

Let us note, moreover, that the inequalities H29|l and 
(|34ll can not be applied to the model of the LMC disk and 
bar we assumed in this paper, because the first property, 
that of homogeneity, is lacking. 

It is therefore useful to derive a more general inequal- 
ity, applicable to all kinds of density distribution. 

6.1. A more general inequality 

We start from equation (|25|l . but we relax the hypothesis 
that the mass density distribution be a function of an 
homogeneous polynomial: 



T < 



47rG' dw psjw) J^^ dz J^^ du piju) 



(35) 



Let us divide in three part the triangular region of 
integration in the plane {u,z), delimited by the bisector of 
the first and third quadrant, by the parallel to the u-axis of 
equation z — w and by the parallel to z-axis of equation 
u = —d: the first region is constituted by the triangle 
delimited by the bisector, the line u = —d and the M~axis; 
the second by the rectangle delimited by the u-axis, the 



z-axis and the lines of equation z = w and u — —d; the 
third by the triangle delimited by the bisector, the z-axis 
and the line of equation z = w. Let us remember that w 
can assume any value between —t and t. In this way, the 
right hand side of Eq. (p?^ is given by the sum of three 
terms: 



47rG 



I -2 J <iwwps(w) 

< /_ddz J^^du pi{u) + /_^dzpi(z)^f^— ^ 



After an integration by part, the first term becomes 
47rG 

dzzps{-z)] 



(36) 



the second term is null, thanks to the symmetry by reflec- 
tion of Ps{z). As regard the third term we observe that for 
any z belonging to the interval (— t, t) it always results 



du p\{u) < j du p\{u) 
^0 



and therefore 

f-t dw P^i'^) Jo d^ lo d^ / .t , , , 

< j_^dwwps[w) 



J ^ dw Ps (w) 



, ■ (37) 

J aw ps 



This way we obtain that also the third term is less than 
zero. In conclusion we obtain the searched inequality for 
the optical depth, valid for any mass density distribution 
of lenses. In particular for lenses in the galactic halo we 
obtain: 



tgh < 



C .10 



dzzpiQii{~z), 



for lenses in the LMC halo we obtain 
AttG /■* 



T'LMCH 



< 



dzzpi,LMCH(^), 



JO 

and for self lensing we obtain 



TSL < 



47rG 



dz z Ps{z). 



(38) 



(39) 



(40) 



Applying this last inequality to the LMC mass density 
distribution used in our (coplanar) model we find 

TSL < 6 X 10-^ 

only about 20% higher than the calculated value. 

Notice that on the right hand side of the inequality H38|l 
Pi.gh(— -z) 7^ Pi.Giiiz), since the galactic halo density is not 
symmetric with respect to the LMC center. The inequality, 
therefore, cannot be obtained by a trivial dropping of the 
factor (1 — z) in the integrand of the expression defining 
the galactic halo optical depth. 
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7. Summary 

The great interest in the location of the observed mi- 
crolensing events towards the LMC is motivated by the 
need to give an answer to the question of their nature. 
Namely, whether (or not) all the events can be attributed 
to known (luminous) populations, so to exclude (or not) 
the possibility for a dark component in the halo in the 
form of MACHOs. 

In this paper we are mainly concerned with the 
possible self-lensing origin of the observed microlensing 
events. In particular we have considered the results of the 
MACHO sur vey. We use the recent pict ure of the LMC 
disk given bv lvan der Marel et al. 1 1|2002|) . and we explore 
different geometries for the bar component, as well as a 
reasonable range for the velocity dispersion for the bar 
population. 

One interesting feature, essentially linked to the as- 
sumed disk geometry, is an evident near-far asymmetry 
of the optical depth for lenses located in the LMC Halo 
(this is not expected, with the possible exception of the 
inner region, for the self-lensing population). Indeed, sim- 
ilarly to the case of M 31 (.Crotts 1992,: Jctzer 1994) . and 
as first pointed out bv lCould I l)l993|) .since the LMC disk 
is inclined, the optical depth is higher along lines of sight 
passing through larger portions of the LMC halo. We show 
that such a spatial asymmetry, beyond the one expected 
from the observational strategy alone, is indeed present in 
the observed events. With the care suggested by the small 
number of detected events on which this a nalysis is based, 
this can be looked at, as yet observed bv ICouldl l|l993|) . 
as a signature of the presence of an extended halo around 
the LMC. 

In the central region the microlensing signatures are 
strongly dependent on the assumed bar geometry. In par- 
ticular, we have studied the variation (that can be as large 
as 50%) in the self-lensing optical depth due to the differ- 
ent geometry of the bar. However, the available data do 
not allow us to meaningfully explore in more detail this 
aspect. 

As a further step in our analysis, we have studied 
the m icrolensing rate. Keeping in mind lEvans fc Kerins"! 
l)2000|) observation that the timescale distribution of the 
events and their spatial variation across the LMC disk 
offers possibilities of identifying the dominant lens pop- 
ulation, we have carefully characterized the ensemble of 
observed events under the hypothesis that all of them do 
belong to the self-lensing population. Through this anal- 
ysis we have been able to identify a large subset of events 
that can not be accounted as part of this population. The 
introduction of a non coplanar bar component with re- 
spect to the disk turns out to enhance this result. Again, 
the small amount of events at disposal docs not yet allow 
us to draw sharp conclusions, although, the various ar- 
guments mentioned above are all consistent among them 
and converge quite clearly in the direction of excluding 
self lensing as being the major cause for the events. 



Once more observations will be available, as will hope- 
fully be the case with the SuperMacho experiment under 
way llStubbs et al. Il2n02^ . the use of the above outlined 
methods can bring to a definitive answer to the problem 
of the location of the MACHOs and thus also to their 
nature. 
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